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We describe results of electronic Raman-scattering experiments in differently doped single crystals 
of YBa2Cu306+a; and Bi2Sr2(Cai,Yi_a:)Cu2 08. The data in antiferromagnetic insulating samples 
suggest that at least the low-energy part of the spectra originates predominantly from excitations 
of free carriers. We therefore propose an analysis of the data in terms of a memory function 
approach which has been introduced earlier for the current response. Dynamical scattering rates 
T{uj) = 1/t{uj) and mass-enhancement factors 1 4- \{uj) = m*{uj)/m for the carriers are obtained. It 
is found that a strong polarization dependence of the carrier lifetime develops towards low doping. 
In B2g [xy) symmetry which weighs out electronic momenta along the diagonals of the Cu02 planes 
the Raman data compare well to the results obtained from ordinary and optical transport. For 
underdoped materials the dc scattering rates in B\g {x^ —y^) symmetry emphasizing momenta along 
the Cu-O bonds become temperature independent and considerably larger than in B2g symmetry. 
This increasing anisotropy is accompanied by a loss of spectral weight in B2g symmetry in the 
range between the superconducting transition at Tc and a characteristic temperature T* of order 
room temperature which compares well with the pseudogap temperature found in other experiments. 
The energy range affected by the pseudogap is doping and temperature independent. The integrated 
spectral loss is approximately 25% in underdoped samples and becomes much weaker towards higher 
carrier concentration. In underdoped samples, superconductivity related features in the spectra can 
be observed only in B2g symmetry. The peak frequencies scale with Tc. We do not find a direct 
relation between the pseudogap and the superconducting gap. 

PACS numbers: 74.72.-h, 78.30.Er, 78.20.Bh, 74.25.-q 



I. INTRODUCTION 



Cuprate systems show a unique doping-temperature phase diagram, and the normal metallic (NC) and supercon- 
ducting (SC) phases evolve from an antiferromagnetic (AF) insulating one upon increasing carrier concentration. For 
a wide doping range AF fluctuations are observed in the NC and SC states. The normal state especially at low 
doping does not exhibit properties of a simple metal, and two if not three qualitatively different new phases have been 
identifled in the temperature range up to some 500 K. [Q 

One of these phases exists between the superconducting transition Tc and a temperature T* and is characterized 
by a pseudogap which is well pronounced in the underdoped range of the phase diagram and fades away in overdoped 
samples. PJsj-p^ Although the phenomenon has been investigated with various methods the critical doping for its 
disappearance is not clear at the moment. Considering its origin several aspects have been discussed without, however, 
winding up with a generally accepted picture. There are two main directions: (i) The pseudogap A* is a precursor of 
superconductivity. Its characteristic energy is then expected to scale with the magnitude of the superconducting gap 
A, and its momentum dependence should exhibit the same symmetry properties. [ p^|jTl|jr^ -pj|| (ii) The pseudogap is a 
signature of the electronic interactions above Tc but is not directly related to the pairing correlations. Then it can well 
show independent scaling behavior and energy scales. pP,pO[] This means that the relationship between the pseudogap 
and the superconducting gap and its evolution with doping is of particular importance for the understanding of the 
cuprates. 
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The underdoped range of the phase diagram is not only remarkable because of the pseudogap but also for the 
increasing anisotropy of the carrier properties in the normal state. This is consistently observed in k-sensitive 
experiments such as angle-resolved photoemission spectroscopy (ARPES) and electronic Raman scattering (ERS) 
p2| , ^ even if the details of the interpretation are not completely clear yet. In order to understand the anisotropy 
electronic interactions with a structure in momentum space are being considered. AF fluctuations, for instance, could 
lead to the observed characteristics as they are strong at low doping and tend to disappear somewhere in the so-called 
overdoped range. In an oversimplified but intuitive way the momentum dependence emerges directly if one recalls that 
nearest and next-nearest-neighbor spins are ordered in a AF and ferromagnetic way, respectively, in an AF lattice. 
This basic idea and physically more elaborate versions of it have been worked out in some detail. [^,^ At the same 
time the microscopic foundations such as Hubbard and t-J models are studied intensively. [p6|-^9| Nesting properties 
of the Fermi surface can also lead to the observed phenomenology. 1 30 

In the present situation it is not only important to compare theoretical predictions and experimental results but also 
to relate results from different experiments to each other. Starting from a microscopical model the respective response 
functions can be calculated directly. In the case of Raman scattering a modified density-density correlation function 
is to be determined which has been done explicitly only for a few special cases. [^sj^olJsil ] However, the calculations 
become increasingly complicated along with the complexity of the electronic interactions studied. On the other hand, 
a more or less quantitative comparison of the results of different experiments is practically impossible in this way. 
Therefore, it seems worthwhile to derive more commonly used functions from the spectra instead. Such a procedure 
is standard in infrared (IR) or optical spectroscopy. There, apart from the conductivity, particle lifetimes t{uj,T) or 
scattering rates F = l/r and mass enhancement factors 1-|-A(a;, T) — m*{uj, T)/mb with rub the band mass are derived 
from the reflectivity |Q following the relaxation or memory function approach proposed by Gotze and Wolfle [Q. 
We will demonstrate in detail below that the method can be adapted for the analysis of the ERS data. We will apply 
it to our results and utilize it even beyond the simple comparison of results from various experiments. It is important 
to note, however, that the different methods cannot be expected to return identical results for the quantities under 
consideration as they measure different correlation functions. The way a probe couples to a system is described by 
the vertex which can well have a non-trivial dependence on energy, for example due to renormalization and 
influence the response function considerably. IR and Raman are two-particle excitations being (essentially) described 
by current-like and density-like vertices which, in most of the cases, can be expressed to a good approximation in 
terms of the one electron energies in the conduction band e(k) as ja oc 9e(k)/9kQ and 7(k) oc 9^e(k)/9kQ9k^, 
respectively. ARPES and electron tunneling spectra reflect the response of single-particle excitations. Therefore, 
only a qualitative comparison can be achieved. Nevertheless better understanding of the interrelation of the different 
experimental methods is very useful as they are complementary in many senses. They couple, for instance, to different 
excitations and have different energy and momentum resolution. A similar step into the direction we want to pursue 
here for Raman has recently been accomplished for the ARPES data in that the electron self energy S(k, cu) has been 
derived [ p5[ . 

In ERS the energy resolution is fairly good and there is also a limited momentum resolution. In addition, the 
density-like coherence factors in the superconducting state allow to clearly distinguish between the occurrence of a 
gap and of phase coherence among the pairs. This leads to qualitatively new conclusions in the analysis of the SC 
gap and the pseudogap. The momentum resolution has its origin in the k dependence of the Raman vertex 7k which 
depends on the polarizations of the incoming and scattered photons and e'^ and, in the most general case, on the 
full band structure e„(k). As outlined above 7k is often approximated by the second derivative of the conduction 
band contracted with the polarization vectors (effective- mass approximation). In any case, 7k can be decomposed 
into its symmetry components 3>^(k) where L denotes the order of expansion, and where the symmetry is indexed 
by /z = Big, B2gT Aig, and A2g. As the symmetries correspond to certain light polarizations a relation between 
momentum and conflguration space is established which allows to weigh out different parts of the Fermi surface or of 
the Brillouin zone with different polarization combinations. p^-|3^ Using this k resolution of ERS the anisotropy of 
the superconducting gap A(k) has been studied in the cuprates. |36|] 

Normal state anisotropics have already been inferred from ERS results a while ago ||2^ without, however, mapping 
them on k space. The static lifetime as obtained from the Raman spectra was compared to ordinary transport 
More recently, we have investigated normal state anisotropics systematically for different doping levels in Bi- and 
Y-based compounds . For momenta along the principal axes (parallel to the Cu-0 bonds in the planes) as opposed 
to those parallel to the diagonals the quasiparticle lifetimes in the static limit {lo = 0) seem to become extremely 
short in underdoped, still metallic and SC samples. In contrast, there is almost no anisotropy in the overdoped range. 
Independent of doping the ERS data at i?2g symmetry compare well with ordinary transport. |p3|| 

One purpose of the present paper is to extend the study to flnite energies and to derive dynamical scattering rates 
from the Raman data. The range w 7^ must be analyzed with even more care than the dc limit since Raman couples 
to almost all excitations in solids. For a derivation of the carrier dynamics other contributions such as phonons or 
magnons must therefore be identified and eventually be subtracted. Fortunately, selection rules and resonance studies 
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using various energies for excitation allow a relatively safe distinction. Nevertheless, we do not think we are already 
in a position to subtract contributions coming from other than carrier excitations. We rather realize that carrier 
excitations dominate at low energies and will therefore confine ourselves to the range up to 1000 cm~^. The method 
for the determination of dynamical properties will be derived in detail. We show that the influence of energies larger 
than 1000 cm^^ leads only to logarithmically small corrections. The obtained data for scattering rate and mass 
renormalization will be compared to the IR results, in particular we will calculate reflectivities for a few examples. 
The careful study of the dynamical response as a function of temperature allows us also to better understand the 
pseudogap state. As a result we can extrapolate normal state behavior observed at T>T* to the range between T* 
and Tc and obtain a good idea of the full influence of the pseudogap on the spectra as a function of temperature and 
doping. This, finally, leads to a comparison of superconducting and pseudogap states. We find clear indications for 
two different energy scales. 

The paper is organized as follows: In sections II and III details of the samples and of the experiment, respectively, 
will be described. The experimental results are compiled in the next chapter where we focus on raw data only. The 
analysis in terms of the memory function will be summarized in section V.A and described in more detail in the 
Appendix. In section V.B dynamical scattering rates and mass renormalization factors will be derived from the data 
presented in IV using the method described in V.A. The discussion (section VI) is split into 4 subsections. In A 
the derived quantities will be analyzed with respect to their low frequency limits and to their symmetry dependence. 
They we will be compared to results from other methods, specifically optical transport, in B. The pseudogap and its 
relation to the superconducting gap will be discussed in detail in C and D, respectively. We present data for two 
groups of compounds and at least four doping levels for each one. The implications will be summarized in section 
VII. 



II. SAMPLES 

The YBa2Cu306+3; (Y123) crystals were grown in BaZrOs cruciblres j4^,^l[]. BaZrOa has been shown to be 
completely inert and to facilitate the preparation of samples with a purity of better than 99.995%. All crystals were 
postannealed in pure oxygen and quenched. Temperatures and oxygen partial pressures were adjusted according 
to the calibration of Lindemer et al. p2[ . The resulting oxygen concentrations were approximately 6.1, 6.5, 6.93, 
and very close to 7.0 for the samples we call antiferromagnetic (AF) insulating, underdoped, optimally doped and 
overdoped, respectively, in the following. There exists a relationship between the oxygen doping and the position 
of the Raman-active vibration of the apex oxygen along the c-axis |}43| . For the metallic samples we obtained 6.53, 
6.93, and 6.96. Given the uncertainty of the Raman-based determination of approximately ±0.05 the results must be 
considered consistent, and they show independently that the oxygen content in the surface layer studied is close or 
identical to the nominal one. In the AF sample the phonon could not be observed for the polarizations studied. The 
magnetically determined respective Tc values (midpoints) and transition widths (10% to 90%) of the superconducting 
samples were 53.5 K (AT^ = 3K), 91.5 K (AT^ = 0.3K), and 87.0 K (AT^ = l.OK). The Tc of the underdoped sample 
when being illuminated by the laser shifts up by several degrees. For the power and the exposure time used the 
saturation limit applies, and the actual Tc is close to 60 K |44[| . 

The Bi-based samples, Bi2Sr2 (Caa;Yi_a;)Cu2 08 (Bi2212), were prepared in ZrO crucibles. In crystals close to 
optimal doping (without Y or, equivalently, x = 1 and appropriately adjusted oxygen content) the resistively measured 
Tc was generally above 90 K with ATc < 2 K. According to transport and Raman results the sample with Tc = 92 K 
is slightly underdoped. If Ca^+ is replaced by Y^+ holes in the Cu02 planes are filled in and Tc is reduced. The 
samples we used contained 38% Y and 100% Y, respectively. The superconducting sample (38% Y) was well in the 
underdoped range of the phase diagram with a Tc of 57 K (ATc = 5 K). Overdoping was achieved by annealing the 
crystals in oxygen at the appropriate partial pressure resulting in Tc's of 81.6 K (ATc = 0.5 K) and 58 K (ATc = 5 K) 
for the samples we used. It has been shown for several compounds that Tc depends in a unique way on the effective 
dopingp as Tc = T™'*'^(1 — 82.6(p— 0.16)^) |^^. This relation is used to determine the respective p's from the transition 
temperatures. Given the uncertainty in Tc and T™"^ the typical error for the doping is Ap — 0.01. The doping levels 
p for our (SC) samples are then 0.09 (Tc = 57 K), 0.15 (Tc = 92 K), 0.16 (Tc = 90 K), 0.19 (Tc = 81.6 K), and 
0.23 (Tc — 58 K). The Tc of the sample we call optimally doped is smaller than that of the slightly underdoped one 
due to a higher degree of disorder. The increase with time of the maximally obtainable Tc is a generally observed 
phenomenon which can be traced back to an improved crystallografic structure and higher purity. Therefore the 
samples are sometimes classified in "generations" . In this nomenclature the optimally doped sample comes from the 
second generation {T™'^^ — 90 K) while the slightly underdoped one is out of the third generation (T™^'^ ~ 93 K). 

Bi2212 is not stable at the stoichiometric composition and there is always excess Bi in the crystals (typically Bi2.i 
instead of Bi2) [^6p^ which is found predominantly in the Ca (Y) position. In addition, a small amount of Sr may 
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be replaced by Ca or Y. In general, it cannot be expected for Bi2212 that a crystal quality comparable to the one 
of Y123 is obtainable. By partially substituting Bi with Pb the modulation along the crystallographic b-axis can 
be changed or even suppressed completely |Q. In the sample we used the Laue pattern still showed an indication 
of a twofold symmetry but it is reasonable to assume that the distortion is smaller than in Pb-free samples. If the 
doping is changed by varying the oxygen content the ratio of the metal atoms plays a crucial role for the stability, 
and underdoped material may become unstable and disintegrate. 

III. EXPERIMENTAL 

The experiments were performed in pseudo back-scattering geometry with an angle of incidence close to 65° and 
the sample surface perpendicular to the crystallographic c-axis. Due to the large real part of the refractive index 
{n ~ 2) the light propagates almost parallel to the c-axis inside of the sample. For the analysis of the scattered light 
we used a double monochromator with single- channel detection. For constant slit width the resolution depends on 
the wavelength of the scattered light. The slits were set at 550 /im, hence the spectral resolution is approximately 
10, 6, and 3 cm~^ in the violet, the green, and the red range of the spectrum, respectively. All spectra are given 
in units of photon counts per second (cps) and are proportional to the cross section. They are corrected for the 
sensitivity of the instrument but not for the frequency dependence of the optical constants. However, n and k do 
not change by more than a few percent in the region of interest . For excitation an Ar+ laser was used, and for 
most of the experiments the line at 458 nm was selected. In a few cases the lines at 476 and 514 nm were used. 
The maximal power was 4 mW in order to keep the laser- induced heating below 15 K. The actual laser-induced 
heating AT and hence the spot temperature was determined in any individual case via the ratio of the Stokes to the 
anti-Stokes intensity. Energy gain and loss spectra fit very well. Therefore, in addition to a well-defined measuring 
temperature, we have strong experimental support for the fluctuation-dissipation theorem to apply. The beam was 
focused to a spot of approximately 50 x 150 /^m^. The polarizations of the incoming and outgoing photons were 
always parallel to the planes. The coordinate system is locked to the Cu-0 bonds with x = [100], x' = [110], etc. All 
symmetries refer to a tetragonal point group. Big and excitations can be projected out separately with x'y' and 
xy polarizations, respectively. Excitations transforming as Aig cannot be accessed independently. All configurations 
with the incoming and the outgoing photon polarizations being parallel project out the Aig contribution with angle- 
dependent B admixtures which have to be subtracted accordingly. Typically, one uses xx — x'y' and x'x' — xy. As 
described in detail in other publications there exists a relationship between electronic momenta and light polarizations 
through the symmetry properties of the Raman vertex: At Big and i?2g symmetry electrons with momenta along 
the principal axes and the diagonals, respectively, are weighed out predominantly [ ^ . If one is interested to sort 
out contributions which transform like A2g it is necessary to use circularly polarized light in addition to the linear 
polarizations pO| , pT| . 

The samples are mounted on the cold finger of a He flow cryostat with the temperature adjustable between 1.5 
and 330 K. The vacuum the sample surface is exposed to is pumped cryogenically and is better than 10~^ mbar. 
For removing surface layers accumulating at low temperature the sample is heated to 250 K once a week. Then 
the experimental conditions are stable and the spectra are fully reproducible. As a general feature in the cuprates 
the scattering intensity at huj 3> fcsT becomes constant. In Fig. 1 we plot the average intensity between 800 and 
1000 cm~^ (corresponding to 1100-1400 K). A variation within ±7% of the mean value is found. Obviously, there 
is no systematic dependence on temperature. This is confirmed by linear regression analysis which produces the 
horizontal lines shown in the figures; they are not the ordinary average of the data. The statistical scatter of the 
data comes from several sources which cannot be controlled sufficiently including the influence of adsorbed surface 
layers. Since there is no trend it is safe to assume that the overall sensitivity varies accordingly without exhibiting a 
spectral dependence. To correct for those small variations all spectra shown below are adjusted to the mean value by 
a multiplicative constant of order 1. 

IV. RESULTS 

Fig. H shows a compilation of Raman results in Y123 at three doping levels, close to 6.0, 6.5, and 7.0. The pure 
symmetry components are plotted in order to specifically clarify the role of the Aig symmetry. As a consistency check 
the differences xx — x'y' and x'x' — xy can be compared. The observed phonon lines are in complete agreement with 
the published literature p^JS^ . The high intensity ratio of the phonons to the continuum independently proves the 
very high crystal quality. The line at approximately 585 cm~^ which is induced by disorder in the chains 0,E3] is 
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present in all samples, even in the undoped and the fully doped ones. This means that a few oxygen atoms are still 
missing in Y123(7.0) and some are left over in Y123(6.0). 

There are several new features of the continuum which have not yet been shown or discussed in the existing literature. 
Particular attention should be paid to the AF insulating sample (Fig. || (a,d,g)) where all intensities of the contimmm 
are generally smaller than in the metallic crystals; the B2g intensity vanishes almost completely. At present it is not 
clear whether or not the residual 0.5 photon counts per second and mW arc intrinsic or due to surface degradation 
or contamination. A similar intensity is also found in the superconducting state at very low frequency at all doping 
levels (see Fig. |^). The cross section at Big symmetry is higher by a factor of 3, but once again approximately the 
same residual intensity is found in superconducting samples close to optimal doping (see Ref. Q). It seems, however, 
qualitatively clear that at least the i?2g intensity is very close to zero in undoped samples. This is plausible since 
the channel apparently samples the carriers at all doping levels studied |^,^,^ . The intensities we found for the 
undoped sample are also considerably smaller than those observed by Katsufuji et al. for comparable La2-KSrj:Cu04 
(LSCO) [Q. The difference between Y123(6.5) and Y123(7.0) is, to within the experimental error, the same as 
already found earlier by Chen and coworkers | |55[ |. 

To make things a bit more quantitative we will focus on parts of the spectra which are not subject to trivial changes 
due to the thermal Bose function, i.e. on the energy range huj 3> fcsT. As demonstrated in the experimental section 
(Fig. |l|) the intensity between 800 and 1000 cm^-'^ does not depend on temperature below 300 K in a systematical 
fashion. For the study of the doping dependences of the intensities at the three main polarizations we will therefore 
compare the average scattering cross sections measured between 800 and 1000 cm' '1 (Fig. I a). The Big intensity 
decreases continuously upon the reduction of carriers while a kind of kink is found for both Aig and B2g symmetry. 
Moreover, the i?2g intensity is peaked at an oxygen content of 6.5. In this context, it should be mentioned that the 
optical constants for visible light depend only mildly on doping and the internal cross section is expected to 
vary by less than a factor of two when changing the oxygen content from 6.0 to 7.0. Although we did not correct the 
spectra for the change of the optical constants trivial effects such as an increase of the intensity on doping can be 
excluded safely as the variation of the intensities goes reciprocally in the B2g channel on the one hand and the Big 
and the Aig symmetries on the other hand. If the Big and the B2g intensities are added the doping dependence of 
the sum is very similar to that in Aig symmetry. This is noteworthy as in lowest order the variation of the magnitude 
of the Aig vertex on the Fermi surface resembles closely that of the sum of the Big and i?2g vertices. The data of 
Bi2212 are similar to those described for Y123 (Fig. || b). The advantage here is that we have also data for strong 
overdoping. The ratio of the Big to the i?2g intensity has been studied in great detail by Naeini et al. [^6| and is 
found to be in full agreement with the data presented here. 

For a clarification of the nature of the excitations it is instructive to study an extended energy range. B2g spectra of 
differently doped Y123 are plotted in Fig. 4. As we have already seen in Fig. 2 the spectra are almost constant in the 
low-energy part. Between 1000 and 2000 cm~^ the slope changes and the cross section starts to increase substantially. 
The same trend is found for the Big symmetry where scattering from two-magnon excitations is superimposed on the 
linear continuum at all doping levels ||l|,^,Q. In contrast, in Aig symmetry the intensity at high energy transfers 
decreases continuously and extrapolates to zero between 10000 and 12000 cm~^ as already observed a while ago in 
Gd2Cu04 [Q. Spectra for Y123(6.0) at energy transfers between 6000 and 15000 cm~^ show similar trends and, 
in addition, structures in the range between 1.5 and 1.8 eV (12000 and 15000 cm^^) j5^. The complete results we 
obtained at high energy transfers will be published elsewhere. 

As already pointed out by Reznik and coworkers [59|] , there is a relatively strong contribution to the cross section 
which is almost independent of doping. Here we find a well resolvable gap for that excitation of some 1500 cm~^ (200 
meV) for Y123(6.0) and Y123(6.5) and one of the order of 1000 cm'^ (150 meV) for Y123(7.0). Superimposed on this 
so far unidentified continuum (see, however, Ref. ||5l[|) there are contributions from scattering off carriers at energies 
down to zero, from spin excitations between 1000 and 4500 cm~^ and, of course from phonons in the range up to 
1000 cm""'^. The gapped continuum has A2g symmetry as it can be observed with comparable intensity at Big and at 
B2g symmetry. This can be checked with circularly polarized light but will not be the subject of this publication. For 
a thorough understanding of inelastic light scattering in the cuprates the origin of all contributions must be clarified. 
It seems, however, safe to attribute the low energy part of the B2g continuum to scattering from carriers as already 
shown earlier j^,^^ and to neglect an interaction with other contributions such as spin excitations. In the following we 
will therefore assume the i?2g continuum for energies below 1000 cm^^ to consist of a single component. Apparently, 
the situation is more complicated in Big symmetry where no superconductivity-induced features can be found at low 
doping ||,||,||,|6|. 

We turn now to the temperature dependence of the spectra and focus on the B2g response in metallic samples in 
the range between strongly underdoped and slightly overdoped. In this context we have to recall that according to 
the fluctuation-dissipation theorem the Raman cross section and, consequently, the number of inelastically scattered 
photons registered per unit time, 7V^-y(a;,T), is proportional to the imaginary part of the Raman response function 

X77(q - o,w,T), 
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N^^iuj,T) = R^^{1 + nBiL0,T)}x';^{q^ 0,u;,T) , (1) 

where {1 + is the thermal Bose-Einstein function. The response function x-y-yi^) — x'-y-yi'-^) + ix'i'vi^) ior 

inelastic scattering of light from electrons in metals is a generalized density-density correlation function pO| and is 
strongly polarization dependent through the vertex 7. To a very good approximation we can take the limit of zero 
momentum transfer, q = 0. For simplicity we will drop the momentum dependence and the polarization (vertex) 
indices 7 and use a dimensionless version of the correlation function in Eq. Hence, all constant factors (trivial 
and nontrivial such as the density of states at the Fermi energy) and units are plugged into R. 

Fig. H shows the spectral response of Y123 in the normal state after division by the Bose factor according to Eq. (|^). 
For clarity two temperature ranges are plotted separately. Above a characteristic temperature T* the variation of 
the spectra is conventional (Fig. ^ (a-c)) in the sense that no intensity anomalies occur: Spectra with a larger slope 
d'x'' {id) / duj in the limes a; ^ lie above those with smaller slope or, more physically, shorter lifetime r. For T < T* 
spectral weight is lost in the range below 800 cm^^ (d-f). At the same time dx"{'^)/du} continues to increase. 
Therefore, spectra at low temperatures cross those at T > T*. For the observation of the effect at higher doping 
levels the samples have to be very pure and well ordered. In materials studied earlier the anomaly associated with 
the pseudogap could not be identified unambiguously In the present samples the anomaly becomes still very 
weak beyond optimal doping but can be resolved clearly in the data. Consequently, in the almost perfectly ordered 
overdoped sample with an oxygen content close to O7.0 the pseudogap can be seen better than in the optimally doped 
one with oxygen vacancies and clusters pl[ . In the same way two- magnon scattering is seen only in Y123(7.0). The 
results obtained for Bi2212 (Fig. H) are in complete agreement with those for Y123. 

Spectra for the superconducting state are plotted in Figs. ^ and ^. Here the data are shown as a function of 
energy units normalized to the respective transition temperatures in order to demonstrate the scaling of the peak 
maxima with T^. The maxima are at approximately 6 and 7 in units of fc^Tc for Y123 and Bi2212, respectively. 
Somewhat unexpectedly but reproducibly p5| , ^ , no pair-breaking peaks are found at the Aig and Big symmetries 
in the underdoped samples (Fig. ^ and p^ ). Moreover, we do hardly find any temperature dependence of the continua 
at these symmetries in the entire range between 10 and 300 K. However, a couple of phonons gain considerably in 
intensity upon cooling 



V. ANALYSIS 



The purpose of this section is to analyze the electronic Raman spectra of metallic samples at frequencies up to 
approximately 1000 cm^^ in a more complete way than done previously. As discussed in the preceding section we have 
good reasons to assume that the cross section in this energy range comes predominantly from carrier fluctuations. 
Hence we can obtain information about the carrier dynamics. To date only the w = limit has been considered. 
For the study of the dynamical response we adopt a method which has been introduced by Gotze and Wolfle |3^] 
for the current-current correlation function and which has been used subsequently for the analysis of infrared data 
1^,^. In this method the relaxation rate F = 1/r is replaced by —iM where Mioj) = uj\{uj) + iT{uj) is the complex 
memory function. r(cL') and 1 + are the (Raman) relaxation rate and mass enhancement factor, respectively. 
The conditions for the applicability are discussed in the paper by Gotze and Wolfle. ||33|] 

We will give a summary of the relevant expressions here and then analyse the data presented above. A complete 
description of the formalism can be found in the appendix . 

A. Model 

The complex Raman response function x{^) terms of the memory function M{uj) is given by 

uj + M yuu) 

hence the imaginary part describing the spectra (see Eq. (|l|)) reads 

^ ^""^ ^ RTTaRF+FmF ■ ^ ' 

We define the function I{uj,T) which is related to the measured Raman spectra N{uj,T) by 



6 



and obtain expressions for the Raman relaxation rate r(cLi) and the mass enhancement factor 1 + 

[I{uj)Y + [ujK{uj)Y 
loK{lo,T) is the Kramers-Kronig (K-K) transform of I{uj,T) hence 



A-(i^) = --P I di , (7) 



and the normahzing factor R is fixed by a sum-rule, 

R=- [ dujl{uj) . (8) 
Jo 

Since in our case I{i-u) approaches a constant in the dc limit and decays approximately as l/oj at high frequencies 
the integral for K{uj,T) (Eq. (^) converges fast as long as the of the upper cut-off frequency ujc is larger than uj (a 
factor of 3 to 5 is satisfactory according to our tests). R depends on tOc only logarithmically. For the study of carrier 
properties a typical value for Wc is the band width, here 1 eV < TilJc < 2 eV . T has the same unit as uj. 
Several remarks seem appropriate: 

(i) Whether or not the choice of Uc was reasonable comes out immediately when looking at the frequency dependence 
of 1 + \{lo). The mass enhancement should asymptotically approach 1 at large energies, in particular it cannot be 
smaller than unity in a system of non-localized carriers. The strongest influence on 1 + \{uj) originates from the 
constant R which, in turn, is the only quantity with a measurable dependence on lJc- So, there exists a second 
constraint for R. 

(ii) One of the conditions for the approach to be reasonable is the resulting spectral dependence of r(a;) and 
1 + A((jj). If these quantities cannot be described by sufficiently smooth functions the physical meaning of the analysis 
is questionable. There are no restrictions as to the high-energy behavior of r(w) and 1 + A(a;). The only requirement 
is a sufficiently rapid convergence of the integrals (0) and (H). 

(iii) Due to the respective vertices the Raman memory function can be quite different from the one for the conduc- 
tivity as intoduced by Gotze and Wolfle |Q . The respective corrections are discussed in the literature . It should 
be noted that the vertices themselves can well have different frequency dependences. 

^iv) The analogy between Raman and conductivity ends strictly at Tc since the superconducting coherence factors 
are density-like and current-like, respectively. This implies that the carrier lifetime can still be derived from the 
optical conductivity in a two-fluid approximation while it is not accessible in a straightforward way in Raman, 
(v) The complex conductivity a{uj) — <7'{uj) + ia"{uj) can be expressed in terms of M{uj) 

ct(w) = ieouli — \ , (9) 

with 

o 

9 ne 

ojIi = . 10 

Here e and m are the electronic charge and mass, respectively, n is the density of carriers and is the permittivity of 
vacuum. It is well known that a sum rule can be obtained for the real part of the conductivity (see also the Appendix) , 
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a'{uj)dijj = n. (11) 







The conductivity sum rule Eq. (|11|) holds strictly, in particular at the transition to the superconducting state where 
the number of Cooper pairs can be determined from the loss of spectral weight at finite frequencies. At the first 
glance, the sum rule given by Eq. (^) indicates the existence of a similar conservation law for the Raman conductivity. 
However, R in contrast to the number of carriers n is not a conserved quantity, because the Raman response function 
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is given by the generalized density-density correlation function [|60| which, as opposed to the usual density operator, 
does not commute with the interaction part of the Hamiltonian. In the superconducting state a completely new 
scattering channel opens up and R can increase considerably. In the normal state, too, there is no reason for R 
to be, for instance, temperature independent. We will see below, however, that R does not change significantly for 
Tc < T < 330 K. In the limit a; — *■ the Drude expression for the dc resistivity po is recovered from Eq. (^), 

Po - T^ro. (12) 

where Fq is the scattering rate in the static limit. 



B. Application to the Raman data 

We calculate now dynamical scattering rates and mass-enhancement factors using the procedure described in the 
previous section. Before that all the phonons have to be subtracted out. This is a straightforward procedure in most 
of the cases as there is little dependence of the phonons on temperature. Only in Big symmetry in Y123 more care 
must be taken, and additional noise resulted from that manipulation occasionally. In B2g symmetry which we are 
predominantly interested in the contributions from vibrational modes to the spectra is generally small and can be 
dealt with easily. 

In Fig. |ll| we plot dynamical relaxation rates and mass-enhancement factors for differently doped metallic samples 
of Y123 which have been determined from the spectra using Eqs. d) and (|). Below 15-30 cm the response is 
linearly extrapolated to zero as inevitably required by the antisymmetry of x"{^)- The spectra beyond 1500 cm~^ 
were extrapolated by a constant since the observed increase does apparently not come from carrier excitations (see 
Fig. U and section VI), the cutoff frequency was set at Tiujc ~ 10000 cm^^. The influence of the extrapolation and of 
the cut off on the magnitudes of F and 1 + A is of order ten percent. The spectral shape at huj < 1000 cm~^ is not 
resolvably affected for hujc > 3000 cm~^. 

For clarity, the temperature ranges above and below T* are plotted separately. F(w,r) has little structure and is 
essentially increasing linearly with frequency. For T > T* the relaxation rates F(a;, T) differ only by a constant offset, 
which is approximately proportional to temperature, and can be well represented by 

F(a;,r) = ac^ + Fo(T) , (13) 

where a is a constant varying between 0.3 and 1 for the samples studied. Fo(r) is the scattering rate in the static 
limit, and at optimal doping one finds ?iFo(T) = phsT. [3 is of order 2 in satisfactory agreement with ordinary 
transport and earlier Raman results . For underdoped samples the static relaxation rate does not extrapolate to 
zero for zero temperature and Fo(T = 0) > 0. If the temperature is reduced below T* a kink starts to develop in 
the spectra at approximately 400 cm^^ which is more pronounced for the underdoped sample (Fig. |ll| (d)) but also 
clearly visible at higher doping. The relaxation rate for w — > 0, Fq, continues to decrease in a similar way as above 
T*, i.e. linear in temperature. For lo > 500 cm~^ the temperature dependence becomes much weaker. 

The mass enhancement 1 + X{u!,T) has little dependence on temperature above T* and has structure only at low 
energies while approaching a value close to unity above 500 cm~^ (Fig- |ll| (g~i))- This demonstrates that the approach 
and, in particular, the normalization through the sum rule Eq. (^ return qualitatively correct results in the far IR 
range as interactions can only enhance the band mass implying X{uj) > 0. The behavior for hto < 20 cm"-'^ is not 
quite clear since the quality of the data is reduced close to the laser line. It appears that there is a weak (logarithmic) 
divergence at w = 0. Since toX and F are K-K related this is a direct consequence of the linearity of T(lu) down to 
uj — 0. Indeed, X{uj — 0) is finite only if T{uj 0) is "rounded off" or more precisely if dT{uj 0)/duj = 0. This, 
however, cannot be decided on the basis of the existing data. For Tc < T < T* a substantial increase of 1 + X{u!,T) 
is found in the small-energy range. 

The results for the Big symmetry are shown in Fig. |l^. In contrast to the B2g symmetry not only the zero-frequency 
extrapolation value and its temperature dependence but also the frequency dependence change. If we assume that 
T{iLj) still varies approximately linearly with frequency according to Eq. (1^) the parameter a changes from 1 in the 
underdoped sample to 0.7 and 0.35 in the optimally doped and the overdoped one, respectively. In addition to that 
F(tj) shows the tendency to saturate for huj > 500 cm^^ in the overdoped sample (Fig. |l^ c) while still increasing 
linearly in the underdoped and the optimally doped samples (Fig. ^ a and b). In the dc limit the relaxation rates 
Fo(T) decrease with T in the overdoped and in the optimally doped samples. In the underdoped one r(w, T) generally 
increases with decreasing T, and the overall magnitude is larger by more than a factor of two. Due to the phonon- 
subtraction procedure and influence of the laser line the results for the mass enhancement become noisy, in particular 
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at small energies. Above some 200 cm"-'^ they are qualitatively similar to those at B2g symmetry at low and optimal 
doping (Fig. |l2|d and e). In the overdoped sample A seems to be close to zero in the whole frequency range (Fig. |l^f). 
Neither for r(w,T) nor for 1 + A(w,T) characteristic changes can be found around T*. 

The B2g results for differently doped Bi2212 samples are compiled in Fig. They exhibit qualitatively the same 
features as those in Y123. Generally, the structures related to the pseudogap are weaker and found at somewhat 
higher energies at approximately 600 cm^^. On account of the larger number of phonons in Bi2212 the data for 



1 + A(ix', T) are less stable than in Y123 (see Fi 
seen as clearly as in underdoped Y123 (cf. Fig. Ill 



and the qualitative change in the pseudogap state is not 
We suspect that the presence of structures related to 



(g-i/ 
and b_3| 

the pseudogap can be suppressed by imperfections in the sample, and it is very likely that Y-doping of Bi2212 affects 
the Cu02 planes stronger than oxygen deficiencies in Y123. In the superconducting state we will observe consistently 
the pair-breaking structures to be less pronounced in Y-doped Bi2212. 



VI. DISCUSSION 



The physical quantities which can be obtained on the basis of the proposed model allow new insights into properties 
of the cuprates. We will first focus on the static (dc) limit and discuss the relationship to other experiments. Later 
on, we will try to arrive at a better understanding of the pseudogap and the superconducting gap and their respective 
evolution with doping and temperature. 



A. The static limit 

Static relaxation rates have already been derived from the Raman data |2^,^,^,^. In order to obtain absolute 
numbers model assumptions as to the spectral shape of the response were used such as a relaxation rate r(a;, T) = 
y/[auj]^ + [ro(T)]2 where ro(T) = r(w = 0,r) |3|||. The frequency dependence was hmited to F while A was set 
identical zero (see e.g. Refs. [|9|65 ). This is a reasonable guess as we find here 1 + A(w) to be constant and close 
to unity for a fairly large frequency range (see Figs. |ll| and 13). In addition, the influence of A is very small in the 
limit CO = (see Eq. (||)). It is important to note, however, that this form of r(cj,T) returns realistic results only in 
the limit w = while failing for uj ^ 0. For finite co the correct frequency dependence is rather described by Eq. ( p^ ) 
which, in turn, does not at all allow to fit the data with A = 0. This detail shows that the proper treatment leads to 
qualitatively different results. 

A reasonable estimate for ro(r) independent of the details of the spectral shape of x"('^) can be obtained directly 
from an extrapolation of the scattering inten sity to = 0, N{uj^O,T), which is proportional to tq{T) — [To{T)]^^ 
in the same way as dx" {uj,T)/duj (see Ref. [Q). If r(ijj) cx a; at large energies (see Eq. (|l^)) the missing constant 
of proportionality is given by the extrapolation value of x"{^~^oo) in any other case, including constant F, by the 
maximal value of x"(cl;). All these numbers are close to 1/2 at the appropriate frequency uji/2- Then, for fiuJi/2 ^ ksT, 
we have a rough estimate of the constant, according to Eq. (Q) i? = 2iV(ajx/2)7 ^^id the zero-frequency extrapolation 
value of the relaxation rate Fg is given by 

nro(T) = ksT^^^^^t^ . (14) 

For the determination of R we always took the average intensity between 800 and 1000 cm^^, A^iooo: which was 
temperature independent in all experiments (see Fig. nh. These arguments rest on quite general principles such as 
the antisymmetry of the response function and illustrate that the low energy part of the Raman spectra follows rules 
which, even if qualitative, provide an idea of the trends and lead already to useful conclusions independent of the 
origin of the intensity: Whenever the response is collision limited, i.e. v^-q <C F with vp the Fermi velocity and 
q the momentum transfer, the slope of the low-energy part of the spectra scales with the lifetime of the respective 
excitation. 

The results for the static limit are consistent with those obtained in the previous section. The temperature depen- 
dence is reproduced almost identically. For the magnitude discrepancies of 10 or 20% are typical. Figs. |l^ and |l^ show 
the results for %Tq{T) as obtained for Y123 and Bi2212, respectively, from both the extrapolation of the intensity 
according to Eq. (|lj) and the extrapolation of the dynamical scattering rate T {uj^T) to lo = Q as shown in Figs. |ll| 
and 1^. The magnitudes have been adjusted to facilitate the comparison. The agreement with the old data is of 
similar quality, and we have to conclude that the relaxational behavior in the dc limit is a stable result. The most 
striking fact is the difference between the Big and the i?2g symmetry being qualitatively similar in the two compounds: 
Towards low doping the Big scattering rate increases, and its slope changes sign from positive to negative. 
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B. Other experiments 



The properties of the carriers have been investigated by several other methods such as ordinary, Hall, and heat 
transport, optical (IR) conductivity, electron tunneling, and ARPES. Most of these experiments are indicative of a 
strong momentum dependence of the relaxation time Tk(u;,T) without, however, converging into a unified picture. 

As demonstrated in Figs. |lj and ^ the scattering rates found here at B2g symmetry in the limit uj ~ show 
a similar evolution with temperature and doping as ordinary transport. For comparison we plotted dc scattering 
rates as obtained from the resistivities through the Drude formula, Eg. (Ill ), represented as energies hTo. For Y123 
all, for Bi2212 part of the resistivities are taken from the literature ||66| , |67|] as are the plasma frequencies in general 
1^,^. For Y123 a-axis data have been used as the i?2g Raman data are insensitive to chain contributions. The 
discrepancies between Raman and transport are smaller than some 30%, and the qualitative change with doping p, e.g. 
the the zero-temperature extrapolation values ro(p, T— >0) and the slopes dro{p,T)/dT of the relaxation rates are 
well reproduced. Independently, very good agreement was found between ordinary and optical transport. []6|j7|j6^ Up 
to this point we apparently wind up with a consistent picture, and the B2g Raman response probes the relaxational 
behavior in the zero- frequency limit as does IR and transport. 

At finite energies, the Raman results and the IR data for T{u!) agree also rather well up to energies of approx- 
imately 1500 cm~^ [H,0j6|]. As in IR the magnitude of the scattering rates does not vary significantly with doping 
implying that the variation of the dc resistivity with doping is apparently not simply an effect of a change in the 
lifetime. It is rather a combination in that the zero temperature off set is due to the lifetime while the increase of 
the slope originates from the change of the plasma frequency (see Eq. (|l^) and Rcf. |Q ) . The differences between 
the IR and Raman results are more pronounced for 1 -I- A(u;,T). The qualitative change with temperature and the 
order of magnitude are still similar (see e.g. [^), however, the range in which the effective mass depends on energy 
is confined to some 500 cm~^ in our experiment while being very wide in IR [6^ . 

Raman and IR can be directly compared when the dynamic reflectivity is calculated from the light-scattering results 
for F and A via the standard formula. In Fig. 16 we show the results for underdoped and slightly overdoped Y123 
using typical in-plane plasma frequencies lo^i obtained from the reflectivity with light polarized along the a axis. As 
the contribution of the chains is suppressed the reflectivity is smaller than for random polarization but both the 
temperature and the frequency dependences are well reproduced, in particular for the overdoped sample. For the 
underdoped one the influence of the pseudogap below T* ~ 200 K manifesting itself as a kink in the IR reflectivity at 
about 800 cm~^ cannot be resolved in the "Raman reflectivity" while being very pronounced in the far IR |69|. Apart 
from details which certainly deserve further clarification ordinary transport, IR and B2g Raman can be considered 
consistent. 

At high energies {huj> 1500 cm^^) Raman and IR cannot be compared any more which is most strikingly demon- 
strated in a plot of the mass renormalization. 1 -I- A(w) as derived from Raman drops significantly below unity for 
huj > 1500 cm~ after a plateau extending over 500-1000 cm' -1 (Fig. 0). This is a clear indication that the cross 
section cannot come from carriers only as already suspected earlier (see Fig. ^ since negative A's are unphysical in a 
metallic system. It should be emphasized that the magnitude of 1 + A is determined independently by the sum rule 
Eq. (^ in the same way as the relaxation rate F which, in turn, is found to be consistent with other experiments. 
Consequently, an absolutely non-trivial result is shown in Fig. |l^. We arrive at the conclusion that there is a lot of 
evidence for Raman scattering at i?2g symmetry to probe carrier excitations in a range up to 1500 cm~^. Beyond this 
range additional channels contribute to the inelastic cross section which do essentially not vary with doping. 

We want to point out again that IR and Raman cannot be expected to return identical results even if the same 
type of excitations is probed. First, different parts of the Fermi surface are weighed out by the two methods. In a 
tetragonal system, for instance, the conductivity is completely isotropic in the basal plane. Apparently, however, the 
B2g Raman response, at least at small frequencies, and ordinary and optical conductivity emphasize similar parts 
in momentum space, i.e. those along the diagonals. Secondly, the bare response functions are renormalized in very 
different ways (vertex corrections) which depend on the probing energies. The exploration of the latter subject is just 
at the beginning. 

It is obvious that the Big data cannot be linked to ordinary transport. If at all their dc limits are reminiscent of 
the Hall resistivity pxy{T,p) ||7C|-|7^, the transverse heat transport (Righi-Leduc effect) or the c-axis transport 
Jt^ . It has indeed been argued that the characteristics of longitudinal and transversal transport are being determined 
by carriers with momenta along the diagonals and the principal axes, respectively [ [7T| , [72{ . Alternatively, spin-charge 
separation can lead to different relaxation times Then the Big rate would correspond to the spin degrees of 

freedom. Indeed, the Big channel samples magnetic excitations becoming more and more dominant at low doping. 
Consistently, in strongly overdoped Bi2212 with essentially no magnetic correlations present the scattering rates are 
almost identical in the two channels (see Fig. ^5|) while the a-c anisotropy is still very large . So there are probably 
two sources for the anisotropy between the Big and the B2g channel in the Raman data, one reflecting the in-plane 
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anisotropy of electronic properties the other one being related to the a-c anisotropy. The in-plane component can be 
compared to ARPES results. There, the scattering rates derived for momenta close to (tt, 0) increase strongly towards 
low carrier concentration |^^. In conventional theory only electronic states lead to ARPES spectra like that, and 
Raman would then just reflect the same strong in-plane anisotropy of carrier properties. It has been shown recently 
that a strong interaction between electrons and spin fluctuations or scattering of electrons on extended impurities can 
lead to such properties |^5[] . It is not clear at the moment how the anisotropy found by Raman can be related to the 
formation of stripes and the crossover from 2D to ID behavior in the pseudogap phase as proposed by Moshchalkov 
et al. l]?^ . If one applies the same argument as for the chains in Y123 the B2g in contrast to the Big response would 
not be expected to be sensitive to properties of stripes running along the Cu-0 bonds. In this scenario the similarity 
between B2g Raman, IR, and transport can only be explained if the overlap between the stipes is substantial and not 
strictly perpendicular to the bond direction. 

We can summarize this section by stating that transport, far IR and B2g Raman at energies below 1000 cm~^ 
can be described consistently in a simple Drude-like picture as long as we limit ourselves to temperatures above 
T* . At low temperature, T <T*, and low doping p < 0.1 the scenario turns out to be oversimplified: (i) In i?2g 
Raman qualitatively new behavior is observed but the related features are weaker than in IR. (ii) Reproducily for 
Y123, Bi2212, and LSCO |p|, p5| , ^ no signature of the superconducting gap can be observed for at Big symmetry 
in underdoped samples with Tc < 0.7 T™°^. There is no significant change in the spectra either at the pseudogap 
temperature T*. At present this observation cannot be linked to the ARPES results where the pseudogap is clearly 
seen for momenta approximately parallel to the principal axes, equivalent to Big in Raman. 



C. Pseudogap 

For underdoped Y123 and Bi2212 i?2g results in the pseudogap state have already been described elsewhere [^||. 
Here, we will additionally focus on the doping and temperature dependence and propose a new method for the data 
analysis. 

The pseudogap state is characterized by a loss of spectral weight in the frequency range between zero and ap- 
proximately 800 cm^^. This decrease of the scattering cross section is clearly seen in B2g symmetry in underdoped 
still superconducting material when the spectra taken at Tc <T <T* are compared to those at higher temperature. 
However, it would be more physical to relate "normal" and "pseudogap" spectra at the same temperature. Since 
the new state below T* cannot be suppressed we try to construct spectra at the respective temperatures. From the 
temperature dependence at T > T* we have a fairly good idea of how the "normal" evolution of the spectra should 
look like: The mass enhancement does not change significantly and the scattering rate r(LL',r) is essentially linear in 
frequency in the range of interest. With decreasing temperature a constant offset is subtracted which is approximately 
proportional to temperature (see Fig. |l^ (a-c) and Eq. (p^)). The constant offset is completely determined by the 
static limit ro(r). Therefore, we take 1 -I- X{uj,T) and T{uj,T) from above T*, leave A untouched, shift F in a way 
that it matches the experimental value of static limit at the desired temperature Fo(T<r*), and plug the quantities 
into the formula for the normal state response, Eq. (^. In this way we obtain an estimate of how the normal response 
Rx"{i^,T) should look like at temperatures in the pseudogap state. The difference between the extrapolated and the 
measured spectra i?Ax"(aj,T) is then a more realistic approximation of the loss in intensity occuring for T < T*. 
The extrapolation procedure can be checked at high temperatures and is found to be fairly reliable (Fig. |l8| ). 

We can now systematically study the dependence of the pseudogap on doping and on temperature. The results, 
i.e. the reduction of cross section, for Y123 and Bi2212 are plotted in Fig. |l^ and |2^, respectively. For underdoped 
material close to Tc the reduction of the response function as a result of the pseudogap has the same order of magnitude 
as the response itself. The change of the spectra in the new state becomes very small for higher doping levels but can 
still be observed even if the Tc of the sample is already reduced by more than 10% on the overdoped side (Fig. pO[ ). 
In Y123 the effect of the pseudogap does not decrease monotonously towards higher carrier concentration and is 
stronger in the overdoped sample than in the optimally doped one. Since we know from magnetization measurements 
that the optimally doped sample we used here is less ordered than the overdoped one due to oxygen clustering ]6l| ] 
it is straightforward to assume that the pseudogap state may escape detection in even slightly disordered material. 
One of the reasons is the pronounced k dependence of the pseudogap as observed by ARPES since scattering from 
impurities leads always to a mixing of different momenta hence averaging out effects confined to small parts of the 
Fermi surface. In a recent ARPES experiment on Bi2212 the effect of disorder due to irradiation with electrons has 
indeed been studied [^: With increasing disorder first the gap in the superconducting state is filled in then the 
electronic states above Tc are being smeared out. In addition, elastic scattering of carriers on defects leads to an 
independent contribution to the Raman cross section which is not at all affected by the new state. This explains 
consistently that the effect of the pseudogap (i) is weaker in underdoped Bi2212 with Ca partially replaced by Y, (ii) 
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is reduced in optimally doped Y123 with oxygen clustering, and, simply because of the continuous improvement of 
the sample quality, (iii) was not observed earlier in Raman. From these considerations, too, we have to conclude that 
the Bi2212 sample with Tc = 92K is already on the underdoped side of the phase diagram. 

While the maximal amplitude of the spectral change i? Ax" {ujps , T, p) is doping and temperature dependent both the 
center of gravity tkops and energy range E* within which the spectra are subject to changes below T* are apparently 
not. We observe them at approximately 200 and at 800 cm^^, respectively. RAx" {T,p) disappears at a doping 
dependent temperature T*(j>). Due to noise and probably also for small instabilities in the measuring conditions T* 
is not very well defined but seems to approach Tc at a not yet clearly determined doping well beyond the optimal one. 

The Raman data fit reasonably well to the results found with other methods as long as the characteristic temperature 
is concerned (Fig. ^l|). However, it is not yet clear to which extent the results can really be compared. The energy 
scales are different and there is actually no microscopic origin identified yet which the various correlation functions 
derive from consistently. The energy Aq found in ARPES close to (tt, 0) has certainly some relation to the normal- 
state gap as it is also observed above Tc. Supposing that in Raman a two-particle correlation function is measured 
giving structures at twice the energy they are observed by a single-particle probe like ARPES the typical energy scale 
of 500 K or 350 cm~^ is reasonable. However, in contrast to Raman where E* is constant Aq oc T* is observed in 
ARPES |80| . A similar discrepancy is found in comparison to electron tunneling experiments | [Tc| , p^ where the energy 
of the normal-state gap scales with that of the superconducting one. In addition to that, the normal-state gap can 
even be observed at very high doping levels and temperatures |11| . This has partially inspired an explanation in terms 
of a charging effect in the insulating layers between the Cu02 planes On the other hand, in all optical methods 
similar energy and temperature scales are observed [p]j82| although some details are different. The kink in T{lu), for 
instance, is at 400-600 cm"-'^ in Raman while it is at 800 cm~^ in IR where we locate E* . 

Recently, Naeini and coworkers have studied LSCO by Raman at various doping levels. They found a strong 
reduction of the Big cross section at low doping independent of temperature (see also Ref. m5j,|5^) and interpreted 
it by Fermi surface destruction as proposed, for instance, by Engelbrecht et al. or Furukawa et al. [Q. Since no 
temperature dependence was found in this study it is apparent that different subjects are being discussed here and 
there which nevertheless can be related through a common interaction such as spin fluctuations in a way suggested 
by Schmalian and collaborators |8^j8^ . 



D. Superconducting gap and pseudogap 



In the current discussion the relationship between the pseudogap and the superconducting gap plays a prominent 
role (see, e.g., Refs. and references therein). Since the Raman response is described by a density correlation 

function coherence factors of type I |Q are relevant. They give rise to the occurrence of pair- breaking features in 
the spectra ||6^,Q thus clearly indicating the onset of phase coherence in the superconducting state. Consequently, 
the existence of a gap or of pairs without phase coherence on the one hand and of coherent pair formation on the 
other hand leads to fundamentally different structures in the spectra. Hence, Raman can provide qualitatively new 
information by clearly distinguishing between a normal state and a superconducting gap in a way more similar to 
Andreev reflection experiments than to ARPES, IR, or electron tunneling. 

The variation with k and the symmetry of the SC gap in the cuprates have been discussed exhaustively, and a 
prevailing dj.2_y2 contribution seems to be widely accepted not only at optimal doping. The spectral dependence 
of the Raman response is in quantitative agreement with the d-wave hypothesis at least close to optimal doping 
1^,^-38 5^j8^,^2l- Complications arise both in the overdoped and in the underdoped range of the phase diagram. 
At high doping levels the spectra seem to be supportive of a reduced anisotropy [po|,pi[ , however, small admixtures 
of an s-wave component can at least qualitatively explain the shift of the peaks ||53||. More importantly, the low- 
frequency power laws clearly show a dominant (1^2 _y2 component. For low doping the situation is more complicated 
as pair-breaking features cannot be found neither at Aig nor at Big symmetry in samples with a substantially reduced 
Tc (Fig. ^ and [lO|) [pp|, |23| , |55| , p6| , p2 93 1. The pair- breaking peaks rather fade away very close to optimal doping [ pS 
In contrast, superconductivity -induced features can be observed reproducibly in B2g symmetry at all doping leve 
For both compounds, Y123 and Bi2212, as well as for LSCO |5q| they can be shown to scale with Tc 



(Fig. |7| and |). For better visualization absolute and normalized peak frequencies for Bi2212 are plotted in Fig. |25^ 
In units of /csTc the peaks are found at approximately 7. In addition, the response increases linearly over an energy 
decade (in contrast to the convex shape in the normal state) with the slope dx" {oJ,T ~0)/duj scaling with ksTc in the 
same way as the peak frequencies (Fig. H). The energies where normal and SC spectra merge hujo exhibit the same 
scaling behavior as hupeak thus demonstrating the spectral shape to be doping independent as well. We conclude, 
therefore, the gap to be d-wave like at all doping levels. Consequently, if (very) exotic k dependences of the gap are 
not considered, the energy scale relevant for the superconducting state is of the order of 8 fc^Tc at all doping levels 
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studied. 

This is not at all a generally accepted conclusion as many other experiments suggest Aq to scale with at best 
in the overdoped range and to continuously increase in the underdoped part of the phase diagram where Tc decreases 
again |l^j9^ , |98|] . This raises the important question of what is actually measured in the various experiments po| , p9| . 
If low energies, E <^ Aq, are considered such as in measurements of the specific heat or of the magnetic penetration 
depth Aah the gap seems to scale with Tc in the whole doping range while typical high-energy probes, E ~ Aq, such as 
ARPES find a decoupling of Tc and of what is usually identified with Aq . In electron tunneling experiments both types 
of behavior seem to be observable, in that close to E — there is scaling |9|] while the position of the conductance 
peaks increases continously towards low doping |l^,Q . It has been argued that the Raman peak frequencies at Big 
symmetry also continue to increase in the underdoped range ||9C| , |9^ resembling the high values found by tunneling 
or ARPES Jp3| . Although most of the recent studies find the Raman peaks to vanish rapidly below optimal doping 
p3| , ^ 56 9^|the frequency positions and, in particular, the interpretation are still a matter of debate. Part of the 
problems is related to the strong sample dependence which is almost completely absent in i?2g symmetry (cf. Figs. |^ 
andU). At Big symmetry we found differences in the peak frequencies of occasionally 100 cm~^ for samples with the 
same Tc and doping level just as a result of different annealing routes or sample qualities. For this reason we will not 
pursue the discussion of the Big response here. 

At present, agreement between the different methods can be achieved for small energies. Here, tunneling, penetration 
depth, and B2g Raman suggest the existence of nodes being characterized by a gap varying as A.(6ip) = CStp with 
Sifi the distance from the node on the diagonal. In addition, the constant C seems to scale with Tc- In satisfactory 
agreement with the data presented here Andreev reflection experiments find the 2C close QUbTc |2^. For the limited 
resolution there are no reliable ARPES data close to the node. Apparently, the B2g Raman data are consistent with 
the low-energy probes not only in the normal state but also below Te- 
la. B2q symmetry we can compare the energy scales of the pseudogap and of the superconducting gap directly. 
In Fig. sp ectra for two doping levels are shown. In both cases pseudogap data for T ~ Tc have been replotted 
from Fig. |l9|. For the superconducting similarly as for the pseudogap state we show the difference between spectra 
at low temperature and those close to Tc, i?[x"(w, T~0) — Tc)], from Fig. 0. In the underdoped sample this 

difference vanishes for energies higher than Tiluq ~ 500 cm^^ well below E* = 800 cm^^. For Y123(7.0) Tiojq and E* 
coincide. At the same time the superconducting peak moves from 220 to 350 cm^^ while the centers of gravity of the 
pseudogap spectra fiujps remain almost unchanged at 200 — 250 cm"-"^. This may demonstrate that we cannot find a 
scaling behavior of the two gaps. E* is apparently an energy scale which is also observed in IR spectroscopy [0| where 
it seems to be as independent of doping as here. Since E* is approximately constant it is tempting to compare it with 
the magnetic exchange energy J ~ 100 meV being of similar magnitude. There is no energy scale which compares 
sufficently well with the magnitude of the spin gap (Aspi„ ~ 160 cm^^ for underdoped Y123 |7^) to motivate an 
identification. fiUps seems to be close in energy but exhibits a different doping dependence. 



VII. SUMMARY AND CONCLUSIONS 



Electronic Raman spectra for the cuprate systems Y123 and Bi2212 have been presented. We study the dependence 
on polarization, temperature T, and doping p. T varies between 10 and 330 K and < p < 0.23 with p = 0.16 being 
the doping at the maximal Tc. For the analysis of the normal-state data the memory function method is adopted 
for Raman. This enables us to calculate dynamical lifetimes r(a;,T,p) and interaction parameters X{lu,T,p) for the 
carriers. 

In the B2g channel emphasizing electronic momenta along the diagonals of the Cu02 planes we find almost quan- 
titative agreement with ordinary transport and IR results for the dc limit and at finite frequencies, respectively. All 
results are consistent. The relaxation rates F = r""'^ do not vary significantly with p. At Big symmetry a strong 
doping dependence of r(a;,T, p) is observed which has not been detected previously. It is reminiscent of the strong 
in-plane anisotropy found with ARPES or derived from transverse charge and heat transport. However, it shows also 
elements of c-axis properties as the anisotropy between Big and i?2g symmetry is quite different in Y123 and Bi2212. 

Below a characteristic temperature T* spectral weight is lost in the B2g spectra for energies Tioj < E* where 
E* ~ 800 cm"-'^ ~ J with J the magnetic exchange energy. By extrapolating the results of the dynamical response 
at T > T* spectra for T < T* can be sythesized and compared to those obtained experimentally. This enables us 
to estimate the full effect of the pseudogap on the cross section and to determine its dependence on doping and 
temperature. E* and the center frequency of the pseudogap do neither depend on p nor on T in the doping range 
the pseudogap can be resolved. The normal-state gap is clearly observed in the range p < 0.20. Its infiuence on the 
Raman spectra is less pronounced than on refiectivity measurements, a fact which might be related to the different 
averaging in the Fermi surface integrals and to different vertices in Raman and IR. 
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In the superconducting state pair-breaking features are observed in B2g symmetry at all doping levels. The spectra 
are indicative of a gap with dx2_y2 symmetry. The characteristic energy scales with the transition temperature in 
agreement with low-energy probes such as the magnetic penetration depth or the tunneling spectra close to zero bias 
and, similarly as derived from Andreev reflection, is found to be close to Sfc^Tc. 

According to these results it seems more likely that the superconducting gap and the pseudogap are completely 
decoupled and that E* is almost doping independent and close to the exchange energy J in agreement with IR results. 
E* is also the maximal energy up to which normal and SC spectra differ at optimal doping or, more physically, up 
to which SC correlations are effective at the highest Tc- The variation of the "gaps" on the Fermi surface is sketched 
in Fig. ^ for two doping levels. Close to optimal doping (Fig. b) the pseudogap is confined to a small region 
in k space, the superconducting gap is well defined all over the Fermi surface and of the same magnitude as the 
pseudogap. In underdoped material the pseudogap is extended over a large region in momentum space. In addition 
to that, the Fermi surface does apparently not exist any more (indicated by dashes in Fig. [2^ a). As a consequence 
the superconducting gap is well defined only close to the nodes where the Fermi surface exists as well. Here, the 
magnitude of the gap scales with Tc- Probably it is meaningless to speak about a superconducting gap for momenta 
along the principal axes. This energy could rather be interpreted as the maximal binding energy of the pairs Ap [ po| . 

Apparently, there are more than two energy scales in the cuprates (subject to a highly fluctuating nomenclature). 
The most prominent and probably best deflned one is the SC gap which is sometimes also referred to as the coherence 
gap Ac and which exists only for T <Tc. There is increasing support for Ac ~ AksTc. The highest scale is set by an 
energy A* |Q or Ap which is identified with the pseudogap or with the pairing energy. They both decrease 

rapidly with doping and have a similar dependence on p. We speculate therefore that A* = Ap. In IR and Raman a 
third almost doping independent energy E* cf. J \s found which becomes effective for T <T* . Finally, there is a gap 
in the spin excitation spectrum ^spin — T* |^J^,Q. The relationship between the various energies is an important, 
so far open issue. 
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APPENDIX: 



The imaginary part of the Raman spectral function x"^(^,r) is obtained from the rate of inelastically scattered 
photons N-y^{uj,T) via Eq. (|l|) as 

i?77{l + nB(w,T)} 

As in the main part nB{uj, T) = {eyiY>{Tiuj / kgT) — 1)^^ is the Bose-Einstein thermal function, the polarization (Raman 
vertex) indices 7 will be dropped, and a dimensionless version of the Raman spectral function x 'wiU be used. 

Since x is a causal function with the imaginary part describing absorption processes ujx" {^) > holds, and x"('^) 
decays at high energies. In the complex z plane we can therefore define a holomorphic function x{'^) by the integral 

X{z)^- / dCf^ (A2) 



for all nonreal z. Without further consequences the Raman response function x(z) may differ by an additive constant 
from the Raman spectral function x{^) defined by the generalized density-density correlation function. The spectral 
function x" {^) is the analytical continuation of x(^) to the real axis, 

X(a>±i,5)=x'(^)±ix"(^) , (A3) 
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and the real and imaginary parts x'(<^) ^^'^ x"('^)j repectively, are related by the Kramers-Kronig (K-K) transfor- 
mation 



^{oj)='^p Tdi^f^^ (A4) 
" /o ? - '■' 



where p denotes the principal value. We define now the Raman memory function M{z) in an analogous way as it was 
proposed for the conductivity by Gotze and Wolfle js^, 

x(=) = . («) 

z + M [z] 

Though the Raman memory function M{z) is different from the conductivity memory function, its analytical properties 
are the same. M{z) is holomorphic for all non-real z and obeys the symmetry relations M*{z) — M{z*) and 
M{z) — —M{—z). Therefore M{z) can be represented by the spectral function M"(w) as 

where M"{ijj) is the analytical continuation of M{z) to the real axis, 

M{uj ± i5) = M'{uj) ± iM"{uj) . (A7) 

For real frequencies lo the symmetry properties of M{z) imply that M'iuj) is an odd and M"{lj) is an even function 
of Lo hence the K-K relation reads 



M'(.) = ^^p / ^ . (A8) 



Jo C 

With the usual notation M(lu) = u;A(a;) -I- iT^uj) both X{uj) and T{ll!) are even functions of uj. 1 + A(a;) = m* /rri}, with 
77ifc the band mass describes a frequency dependent mass enhancement that measures the strength of the interaction 
of the electrons cither among each other or with other excitations such as phonons. In a system of non-interacting 
or elastically scattered electrons A is identical zero. In the high frequency limit A(a;) generally approaches zero in a 
system of free carriers, r(tt') can be interpreted as the inverse of a frequency dependent (Raman) relaxation time, 
V{ijj) — [r(w)]~^. Of course, A(a') and r(a') depend on temperature as well. and V{ijj) are also related by K-K 



transformations as can be seen immediately from Eq. (A8). Though formally similar the memory function M should 



not be confused with self energy E. While M is a two-particle correlation function S is a single particle property in 
a many-body system. 

At the first glance the memory function M{ijj) can be obtained directly from the measured Raman response x"('-^) 



by calculating the real part x'{^) via the K-K relation in Eq. (A4). However, the experimental spectra (see e.g. 
Fig. ^) do not decay but turn out to be constant or to even increase slightly in the mid-infrared range (up to a 
typical band width of the order of an eV) . This requires the introduction of an upper cutoff frequency uJc in the K-K 



integral Eq. (A4) which turns out to influence the results seriously. In order to overcome this problem we deflne a 



new function x(z) 



Since x(x) = 1 — zx{z), xi^) ^^'^ xi^) have the same analytical behavior. For \z\ oo x{z) decays as l/z, and can 
be continuated to the real axis as xi^ ± i^) = x'(^) Obviously, x'(w) = 1 — ujx'{i^), x"{^) — ~^x"(^)j and 

x"(w) < 0. Therefore x'{^) can be expressed through x"{^) as 



X'(a;) = ^p/lc^f^. (AlO) 







We use now the K-K relation for x Eq- (A4) but express x" by x"j 



X (w) = p / 



Jo C 

Jo \ C - ^ 



P/ (1 + 7^— 2 
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This is equal to 1 — lox'{'-^) hence 



1 - wx'(^) = 



^ 9/ .2 

TT 



With Eq. (AlC) we obtain the sum-rule 



1 . 



(All) 



(A12) 



Whenever x"(<^) decays faster than u!~^ the convergence is guaranteed. However, in any practical analysis a cut-off 
frequency uoc has to be introduced which is typically of the order of the band width. This is particularly important for 
Raman scattering since the cross section does not only originate from carriers which the analysis is tailored for. The 
dependence of the integral on Uc is crucially related to the high-energy variation of the response. It is logarithmical 
if x"{'^) decays exactly as or, equivalently, if x"{.^) approaches a constant value. 

The advantage of the new function x becomes apparent when we return to the original goal, i.e. when we derive 
M{ijj) from the Raman spectra. We start by defining a new function I{uj,T) = —Rx"{uJ,T) which can be expressed 
in terms of the measured Raman intensity iV(w, T) as 



/(w,T) 



Niuj,T) 



Lo{l + nBiLu,T)} 



The sum-rule Eq. (A12) fixes the value of the scaling factor R, 

- / dujl{u;) = R 



(A13) 



(A14) 



Additionally, we define K{uj) = Rx'{i^)/ijJ and express it through the K-K transform of x'{^) given in Eq. ( A10| ), 



P / 



m 



(A15) 



As long as u < 0.2 the convergence is very fast and the dependence on ujc correspondingly weak. By separating 
real and imaginary parts we can now determine both the Raman relaxation rate r(a;) and the mass enhancement 
factor 1 + X{lu), 



' A(w) 



[I{u;)]^ + [u;Kiu;W 



(A16) 
(A17) 



We emphasize that ujc is the only adjustable parameter and that not only the spectral shape but also the magnitude 
of r and 1 -|- A are obtained. Due to the specific choice of I{uj) in Eq. ( A13| ) the spectral shapes of T{uj) and 
1 -|- A(w) become almost completely independent of Uc, and the scaling factor R depends only weakly, in our case 
logarithmically, on the cut off. 

The above formulas are general, and the conductivity (y{uj) = c'(cij) +ia"{uj) can be expressed in terms of x{^) [p3[, 



cr(cj) isaujlixi^) 



(A18) 



where ujpi is the plasma frequency. If we plug that into Eq. ( A12) we obtain the well-known sum rule for the real part 
of the conductivity. 



— TT / au! a [oj) = n 



(A19) 



with the definitions of the symbols used in section V. 

We can define a Raman conductivity a-y-y by Eq. (A18) using the Raman memory function. However, due to 
the differences of the vertices the optical and the Raman conductivity should be considered as different quantities. 
Nevertheless they display similar physical properties. The Raman conductivity in terms of the Raman intensity 
I{uj,T) is given by 
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where i?^^ = R-y-y{ujc,T) in the general case, aryy depends on the polarizations and reflects therefore k dependent 
properties. 
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FIG. 1. Average intensities in the constant part of the spectra (800 to 1000 cm~^) for differently doped Y123 at B2g symmetry 
as a function of temperature. The size of the dots is approximately representing the statistical error. The shaded area is the 
variance of the average. 



FIG. 2. Compilation of Raman results for Y123 at characteristic doping levels. Shown are the pure symmetries as indicated 
in the figure. In order to demonstrate the reliability of the subtraction procedure which is required to obtain pure Aig we show 
both XX — x'y and x' x' — xy. Excepting for the overdoped sample (i) which has the strongest orthorhombic distortion the 
agreement is to within the experimental error. In the overdoped sample although twinned there is preferential orientation. 



FIG. 3. Average intensities in the frequency and temperature independent part of the spectra (800 to 1000 cm~^) for 
differently doped Y123 (a) and Bi2212 (b) as a function of symmetry. In Y123 (a), A\g is apparently the direct sum of B\g + 
B2g- The Bi2212 samples come from different sources and "generations", hence the there is more uncertainty in the data. The 
qualitative behavior such as the maximum in B2g symmetry is reproducible. 



FIG. 4. B2g spectra for differently doped Y123 at large energy transfers. 



FIG. 5. -629 spectra for Y123 at different doping levels as indicated for T > Tc- The upper panels (a-c) show spectra at 
T > T* . In the lower ones (d-f) spectra close to Tc < T* are compared to those at T ~ T*. For clarity the spectra of the 
optimally doped (b,e) and the overdoped (c,f) samples have been multiplied by a factor of 2. The shaded areas indicate the 
anomalous loss of spectral weight. 
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FIG. 6. -B29 spectra for Bi2212 at different doping levels as indicated for T > Tc- The shaded areas indicate the anomalous 
loss of spectral weight. 



FIG. 7. -B29 spectra for superconducting Y123 at different doping levels as indicated. 



FIG. 8. -B29 spectra for superconducting Bi2212 at different doping levels as indicated. 



FIG. 9. Raman spectra for underdoped Y123 for T ~ Tc and T <C Tc at all main polarizations. 



FIG. 10. Raman spectra for underdoped Bi2212 for T ~ Tc and T -C Tc all main polarizations. 



FIG. 11. Dynamical relaxation rate hr{uj,T) = h/r (a-f) and mass renormalization 1 + X{lj,T) (g-i) for Y123 as derived 
from the B2g spectra via Eqs. ^ and ^. For all temperatures the samples are not superconducting. 



FIG. 12. Dynamical relaxation rate hr{Lj, T) = Ti/t for Y123 as derived from the Big spectra via Eqs. Q and |6| at T > Tc- 

FIG. 13. Dynamical relaxation rate hr{u!,T) — Ti/t (a-f) and mass renormalization 1 + \{lu,T) (g-i) for Bi2212 as derived 
from the _B29 spectra via Eqs. ^ and ^ (see text) at T > Tc. 

FIG. 14. Static (dc) relaxation rates ?iro(T) for Y123 for B^g (a-c) and -B29 (d-f) symmetry at T > Tc- Open symbols 
represent the data derived via Eq. (|l^, full symbols show the zero-frequency extrapolation values of the dynamical relaxation 
rates plotted in Fig. ^ The data are scaled to each other. Typical differences are of the order of 20% or less. The dashed 
lines in (d) and (e) are relaxation rates calculated from the resistivity through Eq. ^ of comparable samples [Q. The plasma 
frequencies are taken from Ref. ]4S|]. 



FIG. 15. Static relaxation rates ?!,ro(r) for Bi2212 for Big and ,629 symmetry at T > Tc using the same symbols as in the 
previous figure. The dashed lines in (a), (b), and (d) are relaxation rates calculated from the resistivity through Eq. ^ (a) 
and (b) have been measured for the respective samples, (d) is taken from Ref. pi7|. 



FIG. 16. Reflectivities as derived from the _B29 dynamical scattering rates V{uj,T) and mass renormalization 1 + \{uj,T) for 
overdoped and underdoped Y123 at various temperatures as indicated. 

FIG. 17. Mass renormalization 1 + \{io) as calculated from the _B29 Raman response of underdoped Y123 in an extended 
energy range. 

FIG. 18. Extrapolated "normal" spectra for the pseudogap state, T < T* . As demonstrated in (a) the extrapolation 
procedure works quite well above T* where the actually measured spectra can be reproduced satisfactorily. In the pseudogap 
state (b) the deviation from the expected "normal" behavior is very pronounced. 

FIG. 19. The pseudogapas a function of temperature and doping for Y123. Shown in the figure is the difference between the 
experimental and the extrapolated response function RAx" {i^,T,p). 

FIG. 20. The pseudogap as a function of temperature and doping for Bi2212. Shown in the figure is the difference between 
the experimental and the extrapolated response function RAy" {iO;T,p). 
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FIG. 21. Comparison of the characteristic temperatures T* as obtained from different experiments as a function of the 
number of carriers per CUO2 plane. 

FIG. 22. Peak hup^ak and merge huo frequencies for B29 symmetry in differently doped Bi2212 samples. The dashed line is 
a parabolic fit to the peak frequencies resembling closely the doping dependence of Tc as described in Ref. lij) . 

FIG. 23. Pseudogap and superconducting gap at symmetry for two characteristic doping levels in Y123. 

FIG. 24. Variation of pseudogap (dark grey) and superconducting gap (light grey) on the Fermi surface. 
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